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Based on statistical approach we described possible formation of spatially inhomo- 
geneous distribution in the system of interacting Fermi particles by long-rage forces, 
and we demonstrated nonperturbative calculation of the partition function in this 
case. It was shown, that particles interacting with an attractive 1/r potential form 
clusters. Cluster is equilibrium structure, if we suppose that average energy of in- 
teraction of two particles is much less than their average kinetic energy kT. The 
analogy between self-gravitation gas and plasma was shown. The dynamics of clus- 
ter formation was considered with help hydrodynamical and statistical approaches, 
and time of relaxation to equilibrium state was found. 
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I. INTRODUCTION. 



The formation of spatially inhomogeneous distribution of interacting particles is a typical 
problem in condensed matter physics. The types of spatial structures and conditions for 
their formation are determined by type of interaction. A cluster is one from forms of spatial 
inhomogeneity. For gas with attraction of Coulomb type between particles (self-acting sys- 
tem) we can't calculate the virial coefficients (we cannot do it for interaction l/r n if n < 3 



271]). Statistical description of the such system developed in Refs 



4 

on the application of the apparatus of quantum theory of field [12j, [l3|, [l4j, [l5|, [ly, [lTj, [23] 
where it was shown that gas with interacting particles is equivalent to two scalar fields corre- 
sponding attraction and repulsion accordingly with an exponential self-action. The partition 



is based 



function can be representation in terms of a functional integral over these auxiliary fields. 
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The extremal conditions for this functional are nonlinear equations. The spatial distribu- 
tion function, which describes the cluster, is soliton solution of the equations. This approach 
gives the possibility to find the spatial distribution of particles, to calculate the cluster's size 
and to determine the temperature of phase transition in to the state under consideration by 
nonperturbative calculation. 

In the papers (], 7], this approach has been applied to investigation of the system with 
long-range interacting Bose particles. It has been shown that spatial distribution of such 
system is inhomogeneous and has view of finite-size cluster. The radius of cluster, conditions 
and dynamics of its formation have been found both for undegenerated states and for Bose- 
Einstein condensation regime. In this article, based on statistical approach 

naaa 
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demonstrate a nonperturbative calculation of the partition function, we solve the system of 
Fermi particles interacting with an attractive 1/r potential, that is we obtain free energy for 
the such system. We obtain expression for equilibrium radius of cluster and we described 
dynamics of the cluster formation. However, we can find exact results in Boltsman limit 
with first quantum corrections. For all temperatures we propose the approach which gives 
the possibility to find radius of the cluster approximately or to evaluate of one. 

So far as the cluster's radius we can understand as size of spatial inhomogeneity of a 
system, then it is present interest to investigate the general properties of spatial distribu- 
tion function of the systems with Coulomb interacting particles - above-mentioned gas of 
gravitating particles and the system of repulsing particles. This can be done in statistical 
and hydrodynamical approaches. Using hydrodynamical approaches we can find the times 
of relaxation of the systems with Coulomb interacting particles. 

II. STATISTICAL APPROACH TO THE SYSTEM OF INTERACTING 

PARTICLES. 

Let's consider an interacting particles' system being in such conditions when, on the one 
hand, wave's thermal length of a particle can be larger than average distance between them, 
so that it is necessary to take into account type of the statistic, but, on the other hand, 
this length is by far smaller than average scattering length, that lets to describe interaction 
classicall y d isregarding dynamical quantum correlations. The Hamiltonian of such a system 
3, 3, 3 is 
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H ( n ) = X £ s n s - \ X W 88 m 8 n 8 i + U ss in s n s , , (1) 

o — ! —I 

s ss' ss' 

where e s is the additive part of particle energy in the state s (for example kinetic energy 
or energy in an external field), W ss ' and U 8S i are the absolute values of the attraction and 
repulsion energies of particles in the states s and s' respectively The macroscopic state 
of the system is determined by occupations numbers n s . The subscript s corresponds to 



variables that descri 



In Refs. in order to investigate thermodynamical properties of the interacting 



De an individual particle state. 
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llj representation has been used for the 



particles' system a Hubbard- Stratonovich 
partition function in the terms of grand partition function 

Zn = ^nf d ^I DLf I ^exp(-5(e,^V)) . (2) 

where S is the functional which we call an effective free energy which is analogous to an 
action in field theory: 

I* s,s 

+<5]Tln (1 - 5f exp (-(3e s + <p 8 ) cos^) + (JV + 1) ln£, (3) 

s 

£ = is activity, \i is chemical potential, (3 = 1/kT is reverse temperature, s and s' run 
all the states of the system, e is kinetic energy, N is number of particles, 5 = +1 for Bose 
particles and —1 for Fermi particles. The two auxiliary fields tp and ip are introduced corre- 
sponding to attraction and repulsion. The partition function (TSJ) is written as a functional 
integral over these fields. W'^M" 1 , are inverse operators of the interaction: uj~) = 5>Li 

° S,S 1 S,S 1 SS bs 6 

where L s > is such an operator for which the interaction's potential is Green function. 

The integral ([2]) is calculated by the method of "saddle-point" [3, Q across the point 
determined by the functional derivatives 4^ = It = as 

J dip dip 

1 Y^u/-i ^e Vs cos^ n 

l V n-i| ^ a sin ip s n 



and the derivative 



dS* ^ Le^ 8 cos t/i 

— = => V — — - 

<9£ V 1 - cosV>. 



0^V 1 ^ ors ■ = iV+l, (6) 

V 1 - 5£,e^cos^. s v ; 
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where £ s = ^e^ 3 . This set of equations (jUl6]) provides a solution of the many-particle problem 
in the sense that it selects the system states whose contributions in the partition function 
are dominant. The third equation is normalization condition for the particle distribution 
function which is determined by the auxiliary fields. It is obvious that, for given statistics, the 
distribution function depends on the interaction nature and intensity. The cluster formation 
corresponds to particle localization within a limited space. In our treatment, the effect is 
reflected in the behavior of the auxiliary fields and chemical potential. 

For the screened Coulomb or Newtonian potential, the inverse operator may be written 

as 

£ --4^ <A - A2) ' (?) 
where q 2 i s the interaction constant, A is Laplace's operator and A" 1 is the screening length 



0,Q : 



15l . Ilq . Il7j ] . The number of realistic interactions, for which the inverse operator can be 



found and general solution of the set (j4]|6|) can be obtained, is limited. That's why we'll 
confine ourselves to long-range both attractive and repulsive Coulomb potential (oc -^) for 
Fermi gas where we'll find the conditions for the cluster formation and cluster parameters. 

In the continuum approximation, the subscript s runs through a continuum of values in 
the system of volume V. When integrating over impulses and coordinates, we bear in mind 
that the unit of cell's volume in the space of individual states is equal to uj = (27ch) 3 . In 
order to avoid unnecessary complications we shall consider particles without spin only. 



III. CLUSTER FORMATION IN THE SYSTEM OF ATTRACTING 

PARTICLES. 

A. The equations for spatial distribution function for system with Coulomb 

attraction. 

Now we shall consider the system of particles interacting by long-range attraction only and 
we shall demonstrate a nonperturbative calculation of the partition function. For Coulomb 
attraction the inverse operator is known to be W~} = A r 5 rr / , where g is the interactional 
constant and A r is the Laplace operator. 

Let's consider the system in the usual thermodynamical limit: number of particles iV — > 
oo and volume V — > oo with fixed N/V. We consider the effective free energy ([3]) in spherical 
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coordinates: V 
energy 



8R 

for Fermi g 



and, neglecting surface contribution, we can write the effective free 



1 ,v££ i ,v , 
2/3 Jo 4vro 2 cj7o i 



2/3 Jo 4vrg 

+AHn£ = / 
Jo 



1 + exp -p 



2m 



dV 



4r„ 



A 3 



h/2 



+ AHn£, 



(8) 



where r q = 2nq 2 j3, m is a particle mass, A 
and / 5/2 (0 



2irrn 



is wave's thermal length of a particle 



_4 

V5F 



^ / °° c/xx 2 In (l + £e~ x2 ) = E^-l)'* 1 ^ is special Fermi function [ig,|2fl|. Sense of 



the auxiliary field p is next - a spatial distribution's function can be expressed by <p as 



p{R) 



d 3 p ^ ex p(-i) e _ 

MR) ' 



(9) 



(2^)3i + eexp (__2i_) e 

where we used the underintegral expression in Eq.([6]). Let's introduce the dimensionless 
quantity r = R/r q instead of R. Then, the effective free energy (jHJ) (in spherical coordinates) 
can be written as 



S = An r dr 



v 



r h/2 



+ iVln£. 
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The saddle-point equation is Lagrange equation for this functional: 



d\ + 2dip + 2r 3 q df 5/2 (&«>) 



0. 



(10) 



(11) 



dr 2 r dr A 3 dip 

This equation selects the system's states (which are described by the field configuration 
<p(R) or spatial distribution function p[(p(R)]) whose contributions in the partition function 
are dominant. Unfortunately, this equation hasn't analytical solution. But the problem 
simplifies in Boltzmann limit. 



B. The solution in Boltzmann limit and dynamics of the cluster formation. 



Let's consider the limit case £ — > corresponding to high temperature and small con- 
centration, using decomposition of the special function f 5 / 2 (£e p ) in row on orders of the 
activity £. Then the effective free energy([8]) is reduced to the simpler expression 



S = 4ir 



■(Vipf 



A 3 ' 



r 2 dr + iVlnf. 



(12) 
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That is the system with Coulomb attraction is equivalent to single scalar field p(r) with 



231 ]. however therm 



an exponential self-action. Analogous expression was obtained in Ref. 
iVln£, which fixes number of particles, is absent there. The saddle-point equation for the 
functional ( 1T21) is 

rPtn 2r 3 

As it will be shown below, the term can be omitted. In order to connect the auxiliary 

' r or J 

field p with the density p(R) we have to use Eq.Q and to pass to the Bolzmann limit £ — > 0: 

where we introduce the new variable a = exp ((p/2) and let's mark in a 2 = r 3 /A 3 . Then we 
can rewrite Eq. liTBl for the field p as equation for the density as 

d 2 a 1 {® a \ 2 _|_ f a 2 a 3 _ q n§\ 
dr 2 a \ dr J 

This equation has a soliton solution [3j 

A 1 (16) 
y/£a cosh Ar ' 

where A is an integration constant. Any soliton solution corresponds to a spatially inhomo- 
geneous distribution of particles - a finite-size cluster. The corresponding asymptotics are 
cr 2 = 1 for r = d, where d is the cluster size, and a — ► as r — > oo. This solution describes 
the presence of particles in the inhomogeneous formation of the size d and the absence of 
particles at infinity (Fig{TJ). This spatial distribution is compressing to line a = 1 or (p = 
when T — > oo, y — > 0, that's why the field <p — corresponds to spatially homogeneous 
distribution in the system. 

Let's substitute solution ([TBI) in the effective free energy f[T2"j) : 

.2 2 



S = 4:71 J (A — 2£a a j rdr + In £ . (17) 

Then we'll integrate using the decomposition 1/ coshx ~ 1 — x 2 /2 in power series of x = 
Arf « 1: 

S = - V 4^ + Nln Z- ( 18 ) 

The integration constant A 2 is being found from the asymptotics cr 2 (d) = 1: A 2 w £a 2 + 
£ 2 <i 2 a 4 . Thus, we have the result: 

S = -^£ + Ann£--^ 2 dV. (19) 
A 3 A 3 
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FIG. 1: The spatial distribution function cr(r) 2 in a cluster at different temperatures (schemat- 
ically). The solid line corresponds to a lower temperature T\, the dash line corresponds to a 
higher temperature T2. The dot line a 2 = 1 corresponds to mean fixed density of the system and 
determines equilibrium radii of cluster cfoi and ^02 under the above-mentioned thermodynamical 
conditions. 

Assuming that the average energy of interaction of two particles is less than the average 

r 3 N 

kinetic energy ~ kT of a particle -y- < 1, we can find the activity £ from the equation 

X 3 N 2d 2 a 2 \ 6 N 2 



T7T cLS 

<9£ 



£0 + 6,, 16,1 « 6 



> 



(20) 



where £0 and £ 9 are activities of ideal gas and first correction for the interaction accordingly. 
Then integrating Eq.(l2"l) with the effective free energy (1191) on saddle-point (1201) we obtain 
the partition function as 



Z N x exp 



Je ? -ivin (l 



+ — £ 2 d 2 a 2 
£0 J + ^ 



(21) 



where Z Q N is the partition function for ideal gas. Knowing it, we can find a decomposition 
of free energy on cluster's radius d which plays a part correlation length 



Fo-kT 
F + kT 



+ — £ d a 



a 2 \ 3 N' 



d 2 + 



a 



4 X 6 N 3 



-d A 



(22) 



V ™ ' V 2 

where F is free energy for ideal gas. We can see that the d 2 term is always negative 
_ a 2 \^N 2 < q ? faofc j s cluster exists in the system with Coulomb attraction at any thermo- 
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dynamical conditions. Minimizing (|22|) by the size of the cluster 

d 2 a 2 X 3 N 



dF lrT1 da 2 X 3 2N 2 

— - = —kl 

dd V 

we obtain the optimum radius of the cluster, 



1-4- 



V 



0, (23) 



d 2 = (24) 
ANr 3 v ; 



or, for the dimension value D = d ■ r q , 



D ° = 42^N- (25) 

This expression means that equilibrium size of the cluster is defined by balance of the 
two forces. First of them is attraction Coulomb energy, which aspires to compress the gas. 
It is represented by the multiplier q 2 (in sense of interaction constant). The decrease of the 
cluster's size with the increase of average density in the system N/V is connected with a 
closer packing of the particles in the cluster on account of the increase of the energy of inter- 
action. The second of them is the thermal energy, which creates positive pressure resisting 
to the compression. It is represented by the multiplier fcT.With increase of temperature the 
size of the cluster d is rising, however the central density is decreasing. When T — > oo the 
radius of the cluster d — > oo and the central density cr(0) — > 1. Such behavior the value d 
tells us that d plays a part of correlation length which in the point of phase transition (in 
our case is T — >• oo) in infinitely large. Such situation is realized due to the long-range type 
(~ 1/r) of Coulomb interaction. 

On other hand Eq. (125ft may be understood as such such distance in the system on which 
the essential deflection from the average fixed density N/V is observed. If we assuming that 
q 2 = Gm 2 , where m is particle mass and G is gravitational constant, this expression may be 



understood as Jeans length [221 ] . 



Let's consider the dynamics of cluster's formation. For it we'll use the equation of motion 

as 

dd dF 

m = -*W (26) 

where x is a some constant. Applying (|22|) and (|24"I) we have: 

d + rjd 3 - r]d 2 d = 0, (27) 
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where we marked in rj = x 2rf4 T . The solution of Eq. (j2"Tj) is 

2 d l 
d = l + Cexp(-2r7d§t)' (28) 

where C is integration constant. Let's consider the next initial states. 

• Let the radius of a spatial inhomogeneity was larger than the equilibrium size d(t) > d 
in some moment of time t > 0. Then the constant — 1 < C < 0, that means the 
cluster's increase from spatially homogeneous distribution (where d — > oo) to the 
equilibrium size ofo- 

• Let the radius of a spatial inhomogeneity was smaller than the equilibrium size d(t) < 
d . Then the constant C > 0, that that means the cluster's evaporation from initial 
distribution with central density which is larger then equilibrium density at the given 
temperature. 



We can see that size of the cluster approaches to the equilibrium size (1231) asymptotically 
d = do[l+C exp(— 2?7<ioO] when deflections are small. The term \j1r\d\ can be understood as 
time of relaxation. If some fluctuation of density has appeared in the system, then it brings 
about the appearance of the potential's gradient. In its turn, it brings about the spatial 
inhomogeneity - cluster with the size approaching to the equilibrium value asymptotically. 
In others words, cluster is equilibrium structures if we are assuming that the average energy 
of interaction of two particles is less than the average kinetic energy ~ kT of a particle 
-y- <C 1, when the term in Eq. (fl3l) can be omitted, as it will be in the next subsection. 

Since the number of particles in the system iV —>■ oo, but the number of the particles in 
the cluster is finite, then this means that our system disintegrates to the infinity multitude 
of clusters of the size D each. After that the process of cluster formation repeats again, 
where the early formed clusters plays a part of particles. In others words, the free anergy 
of such system hasn't absolute minimum and each state of the system is analogous to false 



vacuum in theory of field 121 ]. 



1. Single-particle motion. 

Let's consider single-particle motion in order to ground our assumptions in the previous 
section for calculation of cluster's size. For evaluation of cluster's size we can make use the 



10 



another method without calculation of the free energy. Let's we assume that the average 
energy of Coulomb interaction of two particles is less than the average kinetic energy ~ kT 
of a particle, when -y- <C 1 that has been named by "local ideality approach" [24J. This 
must lead to the two consequences: 



the activity must be as 



V/N 



(29) 



where £o — is activity for ideal gas. That is the system's activity differences from 
activity for ideal gas little. 



The center density (Fig{TJ) differences from the average density little: p(0) ~ mN/V. 
Then we can think that A 2 ~ ^r?? and we can evaluate the term or -jf in 

VI r or r or 

Eqs.Q13115| with help the solution jTBJ): 



I 2 So- I 
max - -7t— 

I r or I 



dr 2 " or 1 ^ \J 

approach. 



r*N 



<C 1. That is this term can be omitted in the local ideality 



After that let's calculate potential energy of a particle U(r) in the cluster using Boltzmann 
distribution 

U(r) 



a(r) = cr(0) exp 



Then using (fl6l) we have that 



U(r) 



-kT In 



a(0) 



kT 



2kT In [cosh (Ar)) 



(30) 



(31) 



The field U(r) is self-coordinated field of the rest particles in the system (in the cluster). 
Using the assumption (1291) about interactional and thermal energies and considering small 
distances from the center of the cluster as r/r q < 1, we can write that 



„ _ kTNrl _ 
U{r) ss kTA r ss - " r 2 . 



(32) 



We can see that a particle in the cluster is in potential well created by all cluster's particles. 
Then cluster's radius d can be understood as such distant from center of the cluster that 



U{d) = -kT, 



(33) 



11 




FIG. 2: The potential energy U(r) of a particle in the cluster (schematically). The solid line 
corresponds to a lower temperature Ti, the dash line corresponds to a higher temperature Ti. The 
particle does finite motion in the potential well created by the all particles in the system. Then 
the radius of the cluster d we determine as distance from the cluster's center to the point of stop 
of a particle with average kinetic energy \kT. 

where, the multiplier 1/2 corresponds to radial degree of freedom. This means that d is 
distance from the cluster's center to the point of stop of a particle with average kinetic 
energy |&T (FigfSJ). That is radius of the cluster we calculate by self-coordinated manner. 

Then solving the equation ( 1331) we has the equilibrium radius of the cluster as d 2 = ^r^, 
which coincides with Eq. (1241) with exactness to constant multiplier. Since two particles 
attracting by low 1/r, they can not fall on each other in Newton mechanics. Hence spatial 
distribution is regular everywhere in the limits of the local ideality approach, however it is 



not so in the total case 25] . Such mechanism of cluster formation we shall call by the thermal 



mechanism. It is necessary to notice that single-particle motions in cluster is superdiffusive 



26]. 



C. The cluster's size at all temperatures. 

In order to calculate radius of a cluster at any temperatures we must solve Eq. fllll) and 
obtain spatial distribution function of type (fT6j) . We can't do this by cause mathematical 
difficulties, however, if we suppose that gas is ideal locally, that is <C 1, and considered 
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small distances from the coordinates' center R/r q < 1 where \ip\ < 1 (see Eq. (jT4l) and 
Eq. (|16p ). then our problem simplifies. 

For example, let's consider Eq. (TT3"]) and decompose in a row: 

cPm 2r 3 

+ ^ + O(<p) = 0, (34) 
where we suppose that £ « £ - Its solution and corresponding density is 



^--eo^r ^p = -exp I- — 



9 r 2 



Then we can obtain potential energy of a particle in the cluster as 



U(r) = -JfeTln 



exp 



tor, 



A 3 



<? r 2 



kl r, 



(35) 



(36) 



which has be obtained early (1321) . Then, using the condition C/ (cZ) = |/cT and expression for 



the Boltzmann activity £o 



\ 3 N 
V 



we obtain result 



kT 



kid 2 = -kT => d 2 
A 2 



V 



2r 3 iV' 



(37) 



which coincides with Eq.( 

Now let's to calculate the radius of cluster for any temperatures using Eq. tfTT!) and de- 
composing one in a row as 



<9r 



A 3 V <9y? 



(38) 



y?=0 



Analogously, the size of the spatial inhomogeneity is determined by the term at zero degree 
of the field (p. However, the calculations with spatial distribution function (jUj) is somewhat 
problematical because the underintegral expression can not be integrated in elementary 
function. That's why we have to find another method. If we suppose that gas is ideal 
locally, that we have to suppose that motion of particles is quasiclassical [2lJ. Than we 
propose to use the particle's energy in Boltzmann view 

p{r) 



U(r) = -kT In 



P(0) 



-fcTlne^ (r) = -kT<p(i 



(39) 



where the distribution ( 1T41) has been used and the field ip is solution of the Eq.(l38l). however 
the radius of the cluster d we determine as 



U(d) 



(40) 
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where Sp is Fermi energy which is determined by the fixed average density of the system 
AMD - The term \kT is a particle's kinetic energy caused by thermal motion and ^Ep is 
average particle's kinetic energy in degenerated Fermi gas (because all particles can not be 
in the level e = by force of Pauli principle). The term |ep can be understood as repulsive 
"statistical potential" too. Then the assumption about local ideality of gas is written 

2 

as 3? v j N ^ kT + Ep - the average energy of interaction of two particles is less than the 



average kinetic energy of a particle. 

Using the decomposition (|38|) . equation (HUj) and equation (1A3|) for determination of 
the activity £o we can obtain expression for radius of the cluster for Coulomb interacting 
fermi-gas: 



A 3 7V 



(41) 



— - i 3/2 (£o) 

Let's consider some limit cases. In the limit of high temperatures we supposed that 
kT ^> Ep. Then calculation the first quantum correction only in the decomposition (1A3j) we 
can obtain that 

2 _ V 1 A 3 

where \/r q <C 1. For the dimension value: 

2 _ kTV _L h 3 

~ Aitq 2 N + 2 7 /2 7r(27rm) 3 / 2 (fcT) 3 / 2 ' ( ^ 

The sense of the first addendum has been explained in the section flHIBp . The sense of the 
first quantum correction is explained below. 

It is know that so-called " statistical potential" of particles' interaction exists in quantum 
gases jl9| ]. It is repulsion for fermions and it is attraction for bosons. This phenomenon 
is connected with symmetry of particles' wave function. Hence the equilibrium size is de- 
termined by three energies - of Coulomb interaction (the interaction constant is r q or q 2 ), 
thermal kT and above-mentioned "statistical potential" (A or h plays a part of interaction 
constant). That is, in case of Fermi statistics, the quantum energy resists to the compression 
under action of Coulomb force as the small addendum to thermal energy for high temper- 
atures. This explains the sign "+"for the small correction to the cluster's radius. For case 
Bose statistic we can analogous expression but the sign "+" replaced by the sign "— ". 
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In the limit of low temperatures kT Cej? we obtain that 




(44) 



We can see that the equilibrium size is determined by two energies - Coulomb potential and 
the repulsion statistical potential (the constant of interaction is h). That is, the aforesaid 
quantum energy resists to the Colombian compression. The small correction |j means that 
above-mentioned thermal mechanism of cluster formation (as small addendum to Fermi 
statistical repulsion) including if T ^ 0. 

Expression for radius of a cluster at low temperatures can be obtained by the next formal 
means. It is not difficult to notice, that the equilibrium radius of a cluster is determined by 
the multiplier ^p- in the equation (fl3l) after the term at £ ~ £o- m the case T — * we 
can suppose that the radius is determined by above-mentioned multiplier, however because 
temperature is zero, hence kT can not enter into the searching formula. Then it is necessary 
to multiply this coefficient to the value — which cancellates temperature. Then the cluster's 
radius is determined as-^oc-^(j^) = , and we have the formula (jSj). 




IV. FORMATION OF STRUCTURES IN INTERACTING PARTICLE'S 
SYSTEM WITH COULOMB REPULSION 

A. Statistical approach 

Let N + particles repulse by Coulomb low with the interaction constant q 2 . In order to 
ensure stability of such system, we suppose that these particles is dipped in environment 
of contrary sign (with charge of —q), so that the condition of neutralization is executed: 
N + = iV_ = N, but a current created by N + particles only. In other worlds, we consider 
such plasma where constant of interaction is q 2 and charge's carriers of one from the signs is 
immovable. Let's find spatial distribution function for such system, where N — > oo, V — > oo 
and N/V is fixed. 

Let a particle with charge q is situated in the beginning of coordinates. It polarizes the 
gas of particle with the same charge (the environment with contrary sign is immovable). For 
repulsion particles the action has view: 





rv [ (V^) 2 _ J_ 
o [ 4r, A 3 



/ 5 /2(£cosV) dV + N\nZ, 



(45) 
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FIG. 3: The density p(R) of polaron at different temperatures (schematically). The solid line 
corresponds to a lower temperature (Ti), the dash line corresponds to a higher temperature (T2). 
The dot line p = N/V corresponds to average fixed density. 

and corresponding Lagrange equation is 

f± 2r g df 5/2 (£cos?/>) = 

dR 2 A 3 Oib 1 ' 

The term can be omitted as in the case of attracting particles. In Boltzmann limit 
£ — > we can write the action and spatial distribution function as 



/ 



(W>) 2 £ / 
— — cos ip 



Ar q A 3 



dV + N\n£, (47) 



p(it>) = |-cos^, (48) 



with corresponding Lagrange equation 



drib 2r„ . , 

m-^* = - m 

We assume that average energy of interaction of two particles is much less than their 
average kinetic energy, that is £ w £ = ^7^, and we find solution in the area ib{R) G [0, 7r/2] 
where i? G [0, +00] with asymptotics: ij)(0) = 7r/2 and ^(+00) = 0. The first of the 
asymptotics corresponds to absence of particles in the beginning of coordinates, because 
in order to penetrate there a particle must has infinite kinetic energy. Second of them 
means that the concentration equals to average concentration on large distances from the 
center p(R — > 00) = N/V. In other words, the spatial distribution function must be 
pit - such formation called as polaron usually, unlike the system with attraction where 
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spatial distribution is hump - cluster. In other words polaron can be called by anticluster. 
Corresponding solution of Eq. fH9l) is 



tany = tan ^ exp [-J^^-R] • (50) 



4 ~~ 8 ~^ y A 3 
Then spatial distribution function (Figj3]) is written as 



,T*S N 

p{R) = — cos 




7T 



4 arctan | exp | — \j ^ 3 R ] tan — 



(51) 



This distribution is determined by balance of two energies: Coulomb energy which prevents 
from penetration of a particle in the center of coordinates, and the thermal energy (average 
kinetic particle's energy). Then radius of the polaron we can find as distance from its 
center to the point of stop: U(D p ) ~ kT. The function U(r) we can find from Boltzmann 
distribution 

p(^)=p(i^oo)exp(-^P). (52) 

It is using (1511 we have 



U = -kT In 



cos 



4 arctan ( exp ( — \/2!; Q —R ) tan 



(53) 



In order to find radius of the polaron let's apply the equation U(D P ) = ^kT and we obtain 
that 

v « WF < 54) 

This expression for polaron's radius coincides with cluster's radius (|25|) with exactness to the 
constant multiplier. This value we can understand as the size on which the neutralization 
of gas is violated. 

The coincidences sizes D and D p we can explain by the next way. We have seen that the 
spatial distributions in the both cases is result of balance of two energies only: Coulomb 
type and thermal energy. The radii of cluster and polaron are determined by the points of 
stop of a particle in the field of the rest particles. That is mechanism of formation of such 
structures is one and the same - a balance of they two energies. They are differed by roles 
only: in the gas of attracting particles they make spatial distribution of view "hump" and 
in the case of repulsing particles - "pit" Fig.(jl]). 

In the degenerated case situation is analogous, only the Fermi energy plaices a part of the 
thermal energy, that is the equation for the point of stop is U{r) ~ Bp. Let's use the method 
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FIG. 4: The potential energy of a gravitating particle in a cluster - the curves (1) and (1') and 
a particle in a polaron - the curves (2) and (2') at temperatures T\ and T2 accordingly. In the 
Boltzmann limit the radii of the cluster D and the polaron D p is determined as the points of stop 
of a particle with the average kinetic energy kT, which explains these equality. 



described in subsection ( 1III CI) in order to find radius of the polaron. This size is determined 



by the multiplier 4§- in Eq. lHfjj) . but we must get rid of temperature in the final expression. 
For this let's multiply 2r q /X 3 to the dimensionless value kT/ep, where ep is Fermi energy 
( 1A7I) determined by the average density of the system N/V. Then the polaron's radius is 
wrote as 

1 r„ (kT\ q 2 m fN\ 1/3 



This expression coincides with the expression for radius of cluster (jUJ) without the correction 
for temperature in the system of attracting particles. Such coincidence is explained by the 
same as in classical case, but a part of the average kinetic energy kT plays Fermi energy, 
which we can understand as above-mentioned "statistical potential" too. 

B. Hydrodynamical approach 

In order to explain the coincidence of the radii of cluster D and polaron D p we can use 



hydrodynamical approach 22] which lets to look at some general properties of system with 



Coulomb type interaction under other point of view. Let's consider the system of particle 
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which is described be the persistent equation, Euler equation and Poisson equation (" +" for 
gravitating particles and "— " for interaction like electrical charges): 

p— + p(~t = -VP + pV(j) (56) 

A0 = ±4nGp. 



We suppose that perturbations is adiabatical: Pi = fp^i = c 2 pi, where c = ykT/m - the 
speed of sound, Pi ^ P - the small change of pressure. Then the linearized equation of 
motion for fluctuation p is 

( A ~ hw) pi ± k2jpi = °' (57) 



where kj = v/47rGp0 - Jeans wave number. 



c 



Let's described localized oscillation process, when k = 0. Then we have next solution: 
• For the system with Coulomb repulsion the dynamics of increase of the deflection has 



view pi oc exp (iut), where u = ck = J 47r ^ m . That is we have plasma oscillations, 
for self-gravitation system the dynamics of increase of the deflection has view p\ oc 



exp (ut) or p\ oc exp (—ut), where uj = ck = y 47r ^ m .That means exponential increase 
of the deflection or relaxation to the equilibrium state - cluster. 

We can see, that the both process is characterized by one and the same time scale r ~ -/ 
Then we can suppose, that the correspond spatial scale is distant, which a particle cross in 
the process of its thermal motion at the timel/cu: 



LU V 771 



„ <v> kT . kTV 

D « — = = Jj^-^. (58) 



This expression coincides with radiuses D and D p for gravitating particles and repulsing 
particles accordingly. That is both aforesaid systems is characterized by one and the same 
sizes of spatial inhomogeneity. 

As it has been said above, the partition function of self-gravitational system divergences, 
however the states exist where the equilibrium spatial distribution function is solution of 
Eq. (1131) . This means that our system disintegrates (collapses) to infinity multitude of clusters 
of size Dq each, since the number of particles in the system iV — ► oo, but the number of 
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the particles in the cluster is terminal. This multitude is self-gravitating system too and 
it may be collapse. Then the process repeats again. In others words, the free anergy of 
the such system hasn't absolute minimum and each state of the self-gravitating system is 



analogous to false vacuum in theory of field 12]. To evaluate of times of relaxation we can 



using Eq. (158l) . So with each act of collapse, we must substitute m — » nm and iV — > N/n, 
where n is average number of particles in a cluster. Then the time scale 1/ou is constant at 
each process of cluster formation (collapse). 

V. CONCLUSION. 

In this article, based on statistical approach l|, |2j, y, |4| , we demonstrated a nonperturba- 
tive calculation of the partition function, we solved the system of particles with Boltzmann 
and Fermi statistics with long-range interaction of Coulomb type. All results have been 
obtained at the condition that average energy of interaction of two particles is much less 

2 2 

than their average kinetic energy kT or Fermi energy ef- r — *C kT or - i — the 

■ J ^ r yV/N y/V/N 



'local ideality" approach 24]. 



The long-range attraction of particles results to the formation of the cluster of finite size, 
as the initial homogeneous state is instable in the limit iV — > oo, V — > oo with fixed N/V. 
The equation for spatial distribution function has be obtained. However, it is simplified in 
the next limit cases only . 

In case of large temperature (when A <C y/v) the radius of cluster is determined by 
balance of the two energies: the energy of attraction interaction which aspires to compress 
gas, and the thermal energy which creates a positive pressure. This mechanism created 
equilibrium distribution in the the system (fl6l) . in particular, with rising of temperature the 
cluster's size enlarges and its central density decreases accordingly. In the case of degenerated 
system the radius of the cluster is determined by balance of the two energy too: energy of 
interaction which aspires to compress gas, and repulsive Fermi statistical potential creating 
positive pressure. In the total case the equilibrium cluster's size determined by three energy: 
energy of interaction, thermal and repulsive "statistical potential". The size of the cluster 
approaches to the equilibrium size asymptotically in process of its formation. The state of 
the system with spatially inhomogeneous distribution corresponding to the cluster of the 
equilibrium size (]4"3l or 1441) is stable if we are assuming that "local ideality" approach is true. 
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The size of spatial inhomogeneity of two matter: gravitating gas and repulsing particles 
(by Coulomb) is coincides at one and the same interaction's constant in Boltzmann limit 
and degenerated case. This coincidence caused by that, that the spatial distributions in 
both cases is determined by the two energies only: Coulomb type and thermal energy or 
repulsive Fermi "statistical potential". The mechanism of formation of the cluster and the 
polaron is one and the same - a balance of they two energies. This fact can be explained by 
hydrodynamical approach - both matters have same time scale of structure's formation, but 
in the first matter it determines time of relaxation and in the second matter it determines 
of plasma oscillation. This time scales determine spatial scale (the radii of the cluster and 
the polaron) monotonously. 



APPENDIX A: THE EFFECTIVE FREE ENERGIES FOR IDEAL QUANTUM 

GASES. 

Let's obtain some expressions used in this paper for ideal Fermi and Bose gas. Thus we'll 
demonstrate the correctness and greater rationality of our approach as compared with the 
traditional method [191 ]. 

In the case of ideal gas <p = ip = 0. Then the effective free energy ([3]) for the system is 



where 



S= ~^J dv J 47cp 2 dpln{l + £exp{-f3e p )) + AHn£ 



(£) + iVln£, 



(Al) 



h/2 (0 



= / dxx 2 \n(l +te~ x ) 

IT JO V ' 



OO £l 

1=1 1 



(A2) 



is special Fermi function 
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20] ; e p = p 2 /2m is kinetic energy of a particle; A 



ph 2 

2-Km 



IS 



wave's thermal length of a particle. Then the equation for saddle point is as follows: 

1 1 



v A 3 



h/2 (0 



(A3) 



where / 3/2 (f ) = f 



df 5/2(0 „, _ V 



at 



V =N- 
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The partition function (j2j) for this case: 



Z N = exp 



^/ 5 /2 (£)-AHn£ 



(A4) 



where £ is determined from the equation ( 1A3I) . We can find any thermodynamical func- 
tions knowing the partition function. Let's consider the case £ — > corresponding to high 
temperature (the Boltzmann limit), then Eqs. (1A3IIA4|) are reduced to 



Zn = exp 



£- AHn£ 



V N mkT 



N\ \2nh 2 



(A5) 



I 1 

u A 3 



(A6) 



Let's consider the case T — > corresponding to degenerated Fermi gas. Then the special 
Fermi functions are [isj 



/3/2 
/5/2 = 



4 



(ln0 3/2 + -(lnO 



7T 



-1/2 



jj(ln0 5/2 + ^(ln0 1/2 



-or 1 ) 
or 1 ). 



(A7) 



In this case, the expression ( 1A3|) is reduced to 



A 3 
v 



30F 



(InO 



3/2 



where 



6tt 



2\ 3 



(A8) 



(A9) 



2m \ v J 
is Fermi energy. 

We can find analogical expressions for Bose gas. It is necessary to proceed from the 
effective free energy as 



S = ^JdV J 47rp 2 dpln(l - £exp (-/3e p )) + AHn£ 
= --^5/2 (£)+ln(l-£) + Ann£, 



(A10) 



Y^iZi 4j2 - special Bose function 
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20|. It 



where g 5/2 (£) = —fa / °° cfex 2 In (l - £e ; 
is clear that activity always £ < 1 unlike Fermi system; In (1 — £) is effective free energy for 
condensed phase (the addendum with p = is as important as the rest of the sum when 
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Let's find the internal energy U of the system proceeding from the arranged analogy 
between apparatus of thermodynamics in our representation and the theory of field. In 
order to do it let's use the correlation — H = dS mec h/dt and determine conformities H < — > 
U, S mec h * — > Sterm,t < — ► 1/kT, where H is Hamilton's function of the system, S mec h and 
Sterm are the action for mechanics and effective free energy for thermodynamics ([3]) systems 
correspondingly, t is time and 1/kT is reverse temperature. Then, with the help of the 
expression (jSJ), we have 

u -oim = ^^ )= l pv (A11) 

Thus the expressions for thermodynamical function of Fermi and Bose gases obtained by 



our method coincide with ones obtained by usual way 
of the proposed approach. 
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